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We investigate the 2D weakly interacting rotating Bose-Einstein condensate by the tools of quan- 
tum information theory. The critical exponents of the ground state fidelity susceptibility and the 
correlation length of the system are obtained for the sudden change of the ground state when the 
first vortex is formed. This sudden change can also be indicated by the ground state entanglement. 
We also find the single-particle entanglement can be an indicator of the angular momentums for 
some real ground states. The single-particle entanglement of fractional quantum Hall states such as 
Laughlin state and Pfaffian state is also studied. 
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I. INTRODUCTION 

Since the Bose-Einstein condensate (BEG) was ob- 
served in trapped alkali-metal atoms [1, 0, H|, the re- 
sponse of those systems to rotation has attracted con- 
siderable attention. Unlike classical systems, the BEC 
can only gain angular momentum by forming quantized 
vortices when it is stirred. When rotation frequency is 
small enough, no motion of the system can be observed, 
while above some critical frequencies, vortices are formed. 
These vortices are signatures of superfluidity [1] . The for- 
mation of the first vortex is perhaps especially interest- 
ing. Before the formation of the first vortex, the ground 
state of BEC does not rotate while after the formation the 
ground state is a single- vortex state, in which all particles 
rotate around their mass center. Therefore a macroscopic 
symmetry breaking must happen. However, this sudden 
change of the ground state has not been studied quan- 
titatively. Some important physical quantities such as 
critical exponents are still unknown. 

Besides the sudden change of the ground state, many 
authors focused on the the ground state energy and the 
ansatz ground state wave function when the system gains 
a fixed angular momentum through rotation in the weak 
interaction limit [1, i, 0, 0, S M, M, H O- Gener- 
ally speaking, the theoretical methods they used mainly 
include Gross-Pitaevskii (GP) mean field theory and ex- 
act diagonalization (ED). In GP mean field theory, the 
ground state many body wave function is simply ex- 
pressed as a product of N single-particle states, namely 
a non-entangled state. However, sometimes the mean 
field theory cannot give a goo d description to the sys- 
tem. As pointed out in Ref. [14j , when the first vortex is 
formed, the ground state of the rotating BEC will change 
from a product state to a strongly-correlated entangled 
state, leading to the invalidity of the description of the 
mean field theory. Furthermore, when the number of 
vortices N v is large, the mean field theory predicts the 
ground state of the system is a vortex lattice phase [l5[ . 
This prediction is correct only when the filling fraction 
v = N/N v > v c with v c ~ 6 [12j, where N is the particle 
number. When v < v c , the ground state is a strongly cor- 
related vortex liquid phase. Therefore it is the entangle- 



ment between particles that makes GP mean field theory 
invalid. So studying entanglement in rotating BEC will 
be very important to help us understand the properties 
of this system, as what has been done in other condensed 
matter systems [ltjj . However, to our knowledge, this as- 
pect has not yet been investigated. 

In this Article, we use the method of ED in the weakly 
interacting regime to remain the entanglement property 
of rotating BEC. Several tools of quantum information 
theory are used to investigate this system. First we use 
the ground state fidelity and fidelity susceptibility, which 
can precisely locate the critical point of a possibly un- 
known quantum phase transition [l?], [H, QjJ H3, [Hi HI, 
[23j . to study the ground state sudden change when the 
first vortex is formed. By using finite-size scaling analy- 
sis for even TV, we obtain the critical exponents of both 
fidelity susceptibility and correlation length. We find the 
ground state single-particle entanglement can also indi- 
cate this sudden change of the ground state. Then we use 
the von Neumann entropy to calculate the single-particle 
entanglement of the ground states of subspaces of fixed 
z-component angular momentum L z of the system. In- 
terestingly, we find that this single-particle entanglement 
can indicate some L z s of the real ground states, namely 
those stable states. Finally, we study the relation be- 
tween single-particle entanglement and N for some spe- 
cial subspace ground states. We find the single-particle 
entanglement of both bosonic Pfaffian state and bosonic 
Laughlin state diverges logarithmical with N, showing a 
strongly-correlated characteristics of vortex liquid phase, 
while the single-particle entanglement of single-vortex 
state, namely the ground state in the subspace L z = Nh, 
decays with N. 



II. MODEL 

In rotating reference frame, the Hamiltonian of a 
2D rotating iV-boson system with rotation frequency 
£1 trapped in a harmonic oscillator potential is TL = 
YsiLiKoA+U, where H , j = ~^V? + \muj 2 r 2 -VlL z ^ is 
the single-particle Hamiltonian and U oc Yli<j ^( r i ~ r j) 



2 



is the interaction energy. Now we suppose the system 
has a fixed angular momentum L z = LTi{L > 0), then 
its ground state energy when the interaction is absent 
is Eg = (L + N)fko — LMl. The corresponding ground 

state is Y\a—i fi t with a constrain X^i=i h ~ L, where 



= -J= z l e-W I 2 with 



ly is the normalized 



<Pl 

single-particle lowest Landau level(LLL) wave function 
of 7Ycm- ^ the interaction is weak enough (this means in 
Eq.(JTJ NUq < fiw, throughout the calculation we make 
TVlT = hoj/2 to keep the validity of the LLL approxi- 
mation), the dynamics of the system is restricted in the 
LLL, from which other Landau levels are separated by a 
large energy gap 2hw, so that we can use ipi to do the 
second quantization of the Hamiltonian, leading to 



Hl = {L + N)Hlu - LMl + U 



/ Ui ih k,ia\a]akai (1) 

i,j,k,l 



with U i>jtkt l = .J+ 3 jitj\kW. 6i +^ k + l - a ' ( a; ) creates (an- 
nihilates) a particle in the state ipi. The basis Q5l of 
the Hilbert subspace of our system with fixed L z = Lh 
in Fock representation is \NqN\...Nl) with the constrains 
that Yld=a Ni — N and J2f =0 (lNi) = L. Under this basis, 
we can diagonalize the Hamiltonian |T]) to find its sub- 
space ground state l) and the ground state energy 

Now let's consider the total Hamiltonian TL = ® L TLl- 
When fl varies from to w, the one among all I^o.l) with 
the lowest -E"o,i (^) is the real ground state, namely the 
stable state | \&o) 01 "H- The L z of l^o), denoted by L Zt o, 
forms a series of sharp steps from to N(N — 1) (in h 
unit from now on) [3, El- In the first step at SI — 0.75a;, 
L 2j o varies from to TV, corresponding to the formation of 
the first vortex, where the ground state changes suddenly 
from a non-rotating state to a single- vortex state. 



III. GROUND STATE FIDELITY AND 
FIDELITY SUSCEPTIBILITY 

At the beginning we consider a general Hamiltonian 
TL(X) = TL a + XTLb, where A is a parameter that can 
be changed, then the ground state fidelity is defined as 
F = |(* (A + 3A)|*o(A))|. The ground state fidelity sus- 
ceptibility can be calculated from the formula [HI, I2Q, [HJ 



X(A) = - lim 



21nF 



<5A^0 SX 2 



E 



\(^ n (X)\H b \M> (X))f 
[E n (X) - E (X)} 2 



(2) 



where |^ (A)) (|*„(A))) is the ground (excited) state of 
TL(X) and Eq(X) (E n (X)) is the ground (excited) state en- 
ergy. It's obvious that in our system, A = SI and TLb = L z . 
For the system with Hamiltonian TL = @ L Hl, an energy 
level crossing of the ground state exists at SI = 0.75a; 
therefore the ground state fidelity shows a simple drop 
at SI = 0.75a;. One should notice that Eq.Q is valid 
only when there is no degeneracy for the ground state of 
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FIG. 1: (Color online) The lowest three energy levels of TL 
after added by V for N — 5 and TV = 6 (in huj unit). The 
arrow points to the ground state degeneracy when TV — 5. 
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FIG. 2: (Color online) The ground state fidelity F for TV 
5(<5Q = 0.001a;) and N = 4,6, 8(6fl = 0.01a;). 



the system. So to study the ground state fidelity suscep- 
tibility, we have to eliminate this energy level crossing 
first. We adopt the method used in Ref.[14( to add the 
stirring potential V oc — vT) t° t ne Hamiltonian 

TL to generate an energy gap around SI — 0.75a;. Then 
we diagonalize TL which has been added by V in the ba- 
sis 03 = {Jl=q*Bl to calculate x(fi). V must be small 
enough to guarantee that we can still use ipi , the single- 
particle LLL wave function of 7Yo,i in the absence of V 
to do the second quantization. V can be expressed as 

Vo Ei (VW^)aW-2 + V(Z + l)(Z + 2) a\ai + ^j, where 
V < hiu (V = 0.003fio; in this section). We find that 
after adding V, the ground state degeneracy for odd TV is 
not broken completely, meaning the ground state fidelity 
still shows a drop. But for even TV, the ground state is 
non-degenerate so the ground state fidelity is a smooth 
curve (FIG.l and FIG. 2). We focus on even TV in the 
following. 

We calculate x(fi) for even TV. We find that %(f2) 
shows a singularity at SI ss 0.7525a;, where a sudden 
change of the ground state is indicated (FIG. 3(a)). When 
studying the sudden change of the ground state, it's nec- 
essary to consider the problem of scaling. First we will 
study how does the fidelity susceptibility scale with TV 
when SI is at and far from Sl max where the peak lo- 
cates. The height of the peak of x(O) at f2 = O max 
diverges with TV and scales like x(^max) oc TV dc with 
d c w 2.9862 (FIG. 3(b)). When O is far from (either be- 
low or above) Q max , x(O) cx TV d with d = 1 (FIG.3(c)). 
In the thermodynamic limit, Sl max — > f2 c . We are inter- 
ested in the critical exponent /i of the correlation length 
when SI approaches Sl c . It's known that the rescaled fi- 
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FIG. 3: (Color online) The finite-size scaling analysis of x(^) 
for even N. 
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FIG. 4: (Color online) The single-particle entanglement Si of 



the ground state and its derivative 



for N = 6, 10, 14. 



delity susceptibility 
of N»(n - Q max ) [i 



0-x(") 



is a universal function 



FIG. 3(d) shows this function for 
N = 4,6,8,10,12,14 with /i w 1.35. We can see that 
all points of different N locate on a single curve. At 
last we also want to know the critical exponent of x(^) 
when Q approaches £! c . Near f2 c , x(^) must behave like 
oc l/|fi - fi c | Q , where a = (d c - d)/n [H]. So we 
can obtain a w (2.9862 - 1)/1.35 « 1.4713. We also find 
the ground state single-particle entanglement Si, which 
will be defined and discussed in detail in the next section, 
and its first order derivative jMt \ , which shows a sin- 
gularity at Q ~ 0.7525w, can indicate the sudden change 
of the ground state as well (FIG. 4). 



IV. SINGLE-PARTICLE ENTANGLEMENT 

In this section, we want to study the ground state 
single-particle entanglement when the system has a fixed 
angular momentum L z and the stirring potential V is 
canceled. At the beginning, we first introduce the defi- 
nition of entanglement between particles. For a general 
TV-particle pure state |$)i,2,...,at, the entanglement be- 
tween n particles and the other N — n particles can be 
defined as the von Neumann entropy of the reduced den- 
sity matrix of the n particles, namely S n — — Tr(p n In p n ) 



with p n = Tr„ + i ! „ + 2,..., i v(|$)i ! 2,..., J v( < I ) |), where we have 
supposed that |$)i,2,...,jv is invariant under permutation 
of particles (this is satisfied in this paper because of the 
symmetry of boson wave function) so that p n is irrelevant 
with the choice of the N — n particles in the partial trace. 
Therefore to calculate the ground state single-particle en- 
tanglement when the system has a fixed angular momen- 
tum L z , we need to digonalize TLl, numerically under the 
basis *&l z to solve the subspace ground state \^o,l z )- In 
fact |^o,i x ) is expressed in angular momentum occupa- 
tion Fock representation, which should be transformed 
to single-particle state representation. For example, we 
have two identical bosons, one of which has angular mo- 
mentum and the other has angular momentum 1. In 
angular momentum occupation Fock representation, this 
state is |11), but in single-particle state representation, 
this state is -i(i^o)ibl)2+|¥>l)i|¥>o)2) 5 where \(po)x\(px)i 
means the first particle is in the state tpo and the sec- 



ond particle is in the state ipx {tpi 



z l e 



-l-l / 2 ). Af- 



ter completing this transformation, we trace out N — 1 
particles (from the second to iVth for simplicity). Fi- 
nally, the single-particle reduced density operator can be 
expressed in the form pi = Yli=o an d = 

— J2t=o £* m &- O nce L z is fixed, according to Eq.|T|) the 
subspace ground state \^o,l z ) is uniquely determined by 
the operator Y^i,j,k,l^id,k,l a \ a j a hO-l- Therefore Si is ir- 
relevant with Uq. 

We calculate the ground state single-particle entangle- 
ment Si with fixed L z = 1, 2, ...,N(N - 1) + 1. We find 
that Si has a tendency to grow with L z but the whole 
curve shows a behavior of slight oscillation. To see this 
oscillation clearer, in the following we subtract a second 
order polynomial in L z from Si(L z ). For example, for 
N = 8 we subtract -0.0008^ + 0.0817L Z + 0.5057 from 
Si{L z ). In FIG. 5, we find there exist a series of local min- 
ima at L z — L m for Si(L 2 ). Interestingly, most of these 
L m are either the real ground state angular momentum 
L Zt o or the candidate of L z ,q predicted by the composite 
fermion theory Q. Therefore the calculation of Si of the 
subspace ground state \^o^l z ) may give a way of electing 
some L z s of real ground state |^o) (TABLE I). 

Now we focus on the Si of three subspace ground states 
with special L z s. The first one is L z = N, which belongs 
to L Z fi for any N. When L z = N , the ground state wave 
function was conjectured as ipx = JIi=i l( z i~ -z c )e _ ' Zi ' 2 / 2 ], 
where z c — J^iLi z i is the position of the center of 
mass. The single-particle reduced density matrix to ac- 
curacy 0{1/N 2 ) is [| 

r N 

pi(z,z')= / (^Y[dz^j'ip N (z,Z2,...,z N )ip* N (z',Z2,...,z N ) 

J i=2 

The single-particle entanglement can be obtained as 



S^= N w -2/jVln(l/iV) - (1 - 2/jV)ln(l - 2/N). 
FIG. 6(a) shows the S^~ N decays with the growth of 
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TABLE I: In this TABLE, we list L m , where Si(L z ) shows a local minimum; L z ,o, the real ground state angular momentum, 
and Lcf, the candidate of L Zj o predicted by the composite fermion theory Q for N = 5, 6, 7, 8. 
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Lm 


L z ,o 
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6,10,12,15,18,20,24,30 


6,10,12,15,20,24,30 


6,10,12,15,18,20,24,30 
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7,12,14,18,20,22,24,27,30,35,42 


7,12,15,18,24,30,35,42 


7,12,15,18,20,22,24,27,30,35,42 
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FIG. 6: (Color online) (a) Si of the ground state of the sub- 
space of L z = N. (b) red circle: Si of the ground state of the 
subspace of L z = ±N(N - 2) for even N and L z = ±(JV - l) 2 
for odd N, which has a high overlap with bosonic Pfafnan 
state, black cubic: Si of the ground state of the subspace of 
L z = N(N — 1), which is bosonic Laughlin state. 



FIG. 5: (Color online) The single-particle entanglement Si of 
the ground state of the subspaces of fixed L z . A second order 
polynomial of L z has been subtracted from Si(L z ) to make 
the oscillation clearer. The red arrows point to the positions 
of local minima. 



N. This is because when N — > oo, all particles conden- 
sate to the single-particle state (p±. We know that perfect 
condensate has zero particle entanglement [24[ . 

Next we study other two cases which have close rela- 
tion with quantum Hall effect [H, [26|, [2?| ■ One is when 
L z = ±N{N - 2) for even N and L z = |(JV - l) 2 
for odd TV. This subspace ground state has a high 
overlap with the bosonic Pfafnan state of v = 1 Q: 



)pf 



( ) 



\^\'/ 2 . Our nu- 
The other is 



merical result is Si « m(1.1037/V- 1.1369) 
when L z = N(N — 1), which belongs to L z $ for any N. 
This subspace ground state is bosonic Laughlin state of 
v = 1/2: VLau = rfcO* - ^) 2 n^ie- |Zl|2/2 . Our nu- 
merical result is Si w ln(1.941iV-1.081). We can find the 
single-particle entanglement analytically for some quan- 
tum Hall effect states in spherical geometry obtaining 
Si Lau = ln(27V - 1) for v = 1/2 and Sf f = ln(N - 1) for 
v = 1 i 28( . Now we are dealing with a system in disk ge- 
ometry so our results are not exactly the same with this. 
However, the difference is not very large (FIG. 6(b)). 

It's known that when and N v are both large and our 
system has a fixed L z , if v = N/A^ = N 2 /(2L Z ) > v c 



with v c ~ 6, the ground state is a vortex lattice state, 
otherwise the ground state is a vortex liquid state. In 
the vortex liquid regime for v — 1/2 and v = 1, the re- 
sults above show that in the thermodynamic limit Si of 
the ground state is logarithmical divergent with TV. But 
what's the relation between Si of the ground state and 
A^ for some v in the vortex lattice regime? Considering 
in this regime the mean field theory describes our sys- 
tem well, we conjecture that when N — * oo, Si will not 
diverge with A^ for some v in the vortex lattice regime. 

We hope the ground state sudden change when the first 
vortex is formed and the entanglement in it can be inves- 
tigated by experiments. To achieve this goal, an energy 
gap between the ground state and excited states needs to 
be generated by stirring potential (see FIG.l). Consider- 
ing this energy gap decays with the particle number, the 
number of particles in experiments should be restricted. 
On the other hand, for most experiments of rotating BEC 
realized in the laboratory, the filling factor v ~ 10 3 , well 
inside the vortex lattice phase. Experimentalists have 
elaborate techniques to observe vortex lattice, while how 
to detect the vortex liquid state is still a challenge. How- 
ever, some experimental methods have been proposed 
[29l ]. We hope the strongly-correlated characteristics of 
bosonic Laughlin and Pfaffian state as we show here 
by their logarithmical divergent single-particle entangle- 
ment can be verified in experiments in the future. 
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V. SUMMARY 

In summary, we investigate a 2D rotating BEC by tools 
of quantum information theory. The critical exponents 
of ground state fidelity susceptibility and the correlation 
length are obtained for the ground state sudden change 
when the first vortex is formed. We find the single- 
particle entanglement S% of the ground state can be used 
to detect this sudden change. We also find a novel prop- 
erty that S± can indicate some angular momentum L z of 
the real ground states, namely those stable states. At 



last, Si of the ground states in some special subspaces 
of fixed L z are calculated to show the strongly-correlated 
property of vortex liquid phase. Thus some basic proper- 
ties underlying in rotating BEC are clarified and viewed 
from the point of quantum information. On the other 
hand, the formation of the first vortex in rotating BEC 
may provide a signature for the macroscopic entangle- 
ment 30]. 

This work is partially supported by 973 grant No. 
2010CB922904. 
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